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We show that 5-dimensional Kaluza-Klein graviton stresses can slow the decay of shear anisotropy 
on the brane to observable levels, and we use cosmic microwave background anisotropics to place 
limits on the initial anisotropy induced by these stresses. An initial shear to Hubble distortion of 
only ~ lQ~^Q,oh^ at the 5D Planck time would allow the observed large-angle CMB signal to be a 
relic mainly of KK tidal effects. 



Recent developments in string theory have inspired the 
construction of braneworld models, in which standard- 
model fields are confined to our 3-brane universe, while 
gravity propagates in aU the spatial dimensions. A simple 
5D class of such models allows for a non-compact extra 
dimension via a novel mechanism for localization of grav- 
ity around the brane at low energies. This mechanism is 
the warping of the metric by a negative 5D cosmological 
constant . These models have been generalized to ad- 
mit cosmological branes Q , and they provide an interest- 
ing arena in which to impose cosmological tests on extra- 
dimensional generalizations of Einstein's theory |3|-|^] . 

Modifications to general relativity in the direction of 
a quantum gravity theory need to be consistent with in- 
creasingly detailed cosmological observations. The pre- 
mier cosmological test is provided by cosmic microwave 
background (CMB) anisotropics. A detailed calculation 
of CMB anisotropics predicted by braneworld models is 
complicated by the need to solve the full 5D perturba- 
tion equations, involving partial differential equations for 
the Fourier modes. Up to now, only qualitative or spe- 
cial results are known 1^-0] ; in particular, the Sachs- 
Wolfe effect has not yet been calculated because of 5D 
effects Q . The 5D effects are carried by so-called Kaluza- 
Klein (KK) massive modes of the graviton, which can 
generate anisotropy on the brane. In view of the great 
complexity of the full 5D problem, it is worth explor- 
ing partial aspects of the problem. In this spirit, we 
impose some physically reasonable assumptions on the 
5D KK effects in order to estimate the CMB large-angle 
anisotropy. Using COBE observations, this then leads to 
constraints on the KK anisotropy. We find that the CMB 
imposes significant limits on the initial anisotropy. It is 
even possible that the observed large-angle anisotropy 
derives from anisotropic KK gravitational effects. 

Extra-dimensional modifications to Einstein's equa- 
tions on the brane may be consolidated into an effective 
total energy- momentum tensor [^J^ : 

G^iy — K^T^5 = (Tp^ + T^°^ + r^^) . (1) 



Since the brane cosmological constant A is negligible at 
the early times that we consider, we choose the bulk cos- 
mological constant so that A = 0. The local effects of 
the bulk, arising from the brane extrinsic curvature, are 
encoded in T]^^ ~ {T^^)^/X, and are significant at high 
energies above the brane tension Q 

p > A > 10« GeV^ . (2) 

The nonlocal bulk effects, arising from tidal stresses im- 
printed on the brane by the bulk Weyl tensor [^,||, 
are the KK modes carried by Tj^^. We are inter- 
ested here in the astrophysically relevant case of small 
anisotropy and inhomogeneity, and we adopt a long- 
wavelength velocity-dominated approximation, so that 
inhomogeneous scalars vary slowly with position and 
time-derivatives dominate over spatial derivatives. In 
this approximation, we neglect the acceleration and vor- 
ticity of the comoving 4- velocity u'^. Furthermore, since 
we are interested in how anisotropic stresses from bulk 
gravitons affect the shear anisotropy cr^^ on the brane, we 
neglect the matter and KK energy fluxes, and the trace- 
free anisotropic matter stress. The effective total energy 
density, pressure and anisotropic stress are therefore |3[| 

P°*^ = p(1 + p/2A + /Vp) , (3) 
- p (1 + p/X) + p (p/2A + p^^/3p) , (4) 

7rV = 0. (5) 

The brane energy-momentum tensor separately sat- 
isfies the conservation equations, VT^j/ = 0, and the 
Bianchi identities on the brane imply that the effective 
energy-momentum tensor is also conserved: VT°^ = 0. 
In general relativity, anisotropic stresses slow the decay 
of shear anisotropy ||^. Without anisotropic stress, this 
slower decay is impossible in general relativity, but in 
the braneworld, KK graviton stresses can play a similar 
role to anisotropic stress, as we show below. With our 
assumptions, the conservation equations (see |^) reduce 
to 
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p + e(p + p) = 0, 
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where tt^ is transverse as well as tracefree, a dot denotes 
a comoving time derivative, and Q = 3a/ a, where a is an 
average scale factor. 

We also assume that the spatial curvature may be ne- 
glected. This leads via the Gauss-Codazzi equations on 
the brane to a shear propagation equation and a 
Friedmann-likc equation: 
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In Eqs. (|)-(|), there is no evolution equation for 
the nonlocal KK anisotropic stress ttJ^J^. This is the 
anisotropic stress imprinted on the brane by the 5D Weyl 
tensor, and this nonlocal field cannot be predicted by 
brane-bound observers since it includes 5D gravitational 
wave modes. The 5D field equations must be solved to 
determine the brane dynamics completely, and this also 
involves choosing boundary conditions in the bulk. On 
the other hand, starting from a brane-bound viewpoint, 
any choice of ttJ^J^ that is consistent with the brane equa- 
tions, will correspond to a bulk geometry, which can be 
locally determined in principle by numerical integration 
(or approximately, close to the brane, by Taylor expand- 
ing the Lie-derivative bulk equations given in [^). How- 
ever, numerical integration is very complicated (see px| ] 
for the black hole case). Even if it can be successfully per- 
formed, it will not give the global properties of the bulk. 
The bulk geometry that arises for a given ttJ^^ may have 
unphysical boundary conditions or singularities (e.g., the 
bulk corresponding to a Schwarzschild black hole, with 
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■kt/i = 0, has a string-like singularity and a singular 



Cauchy horizon |1 1|). 

We have no exact bulk solutions to guide us in a study 
of cosmological anisotropy. The only relevant known so- 
lution 13 1 is the Schwarzschild- anti dc Sitter bulk that 
contains a (moving) Friedmann brane, which is the ex- 
actly isotropic and homogeneous limit of our case, with 
= 0. In the absence of exact or numerical 5D so- 
lutions, we are forced to make assumptions about the 
KK anisotropic stress tt^ in order to estimate its impact 
on the shear anisotropy of the brane. These assumptions 
should be consistent with the brane equations above, and 
physically reasonable, and are discussed below. 

Observational constraints on the KK stresses arise 
from big bang nucleosynthesis and from COBE measure- 
ments of large-angle CMB temperature anisotropics. The 
KK energy density on the brane introduces a new radia- 
tive degree of freedom at nucleosynthesis. Assuming a 
maximum of 0.3 of a 2-component neutrino species, and 
helium limits of 0.228 to 0.248, this gives (compare p2t^) 
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^ ~s(ii,) < 10-5, s = yi;::^, (n) 



where tjs is the time of last scattering. By Eq. 
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2k^P [1 + p/2A + p'^Vp] /e 
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and then Eq. (|l0| ) gives the nucleosynthesis limit 

S(tns)<0.13. (13) 

For small anisotropy, the volume expansion of the uni- 
verse is determined by the isotropic matter source [we 
assume that p^^ = in the background isotropic so- 
lution of Eq. (H)]. This amounts to treating a to low- 
est order as the scale factor for the isotropic Friedmann 
braneworld. In this approximation, and with equation of 
state p = (7 — l)p, with 7 constant, Eqs. (^ and (||) lead 
to 

(2i + /3)/[7t(t + /3)] , p = A/[3K^jh{t + /3)] , (14) 
sTa/k < 10-9 sec , (15) 

where the bound follows from Eq. (^. This is in agree- 
ment with the Friedmann brane solutions given in [^ . 
(If we take into account the nonlocal energy density in the 
background solution, i.e. p^^ 7^ 0, we can generalize this 
solution when 7 = |: O is still given by Eq. (Q), with /? 
replaced by /3 = /3[1 + 2p^^ / (n'^Xp)], while p ^ (/3/^)p.) 
The low-energy limit is /3 — > 0, when the general relativ- 
ity solutions are regained. 

Cosmological anisotropy on the brane has been con- 
sidered in recent papers. A qualitative description of the 
role of nonlocal anisotropic stress in cosmology is given 
in , and perturbative analysis on large scales is devel- 
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The large-angle CMB anisotropics are constrained by 



], but without finding or assuming a form for 
4IJT5I, Bianchi I dynamics on the brane, with 
vanishing spatial curvature, is studied. In particular, the 
local extra-dimensional modifications to general relativ- 
ity introduce a novel feature to early-time dynamics [p^ : 
instead of shear domination, there is matter domination 
at early times, and the relative shear anisotropy s is a 
maximum when p/A = (2 — 7)7(7 — 1). 

Here we investigate the evolution of shear anisotropy in 
an inhomogeneous universe, at early times (when A may 
be neglected) and in the absence of anisotropic stress 
from fields on the brane, i.e. tt^^ = 0. We need to 
determine the effect of KK anisotropic stress on 
braneworld shear. In principle, the 5D Weyl tensor is de- 
termined by the solution of the 5D field equations, and 
its projection onto the brane then determines ttJJ^. In 
this way, the 5D solution determines an effective evolu- 
tion equation on the brane for tt^. In practice, we know 
of no way to find or estimate this 5D determination of 
TT^. We emphasize that tt^ is a 4D quantity, which is 
an effective 4D anisotropic stress, even though its evolu- 
tion is governed by 5D graviton dynamics. Thus it seems 
reasonable as a first approximation to assume that 
behaves qualitatively like a general 4D anisotropic stress. 
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This should be general enough to cover a wide range of 
bulk graviton effects. Then we can estimate the KK effect 
on large-angle CMB anisotropics. (Qualitative estimates 
in the absence of 5D solutions have also been used to es- 
timate CMB anisotropics in different braneworld models, 
where gravity is modified at large scales rather than high 
energies ]T6|.) 

There is an ansatz that includes all examples of 
anisotropic stresses studied in relativistic cosmologies 
and is physically well motivated. According to this 
ansatz, in the large-scale velocity-dominated approxima- 
tion, the time evolution of tracefree anisotropic stress 
is proportional to the energy density of the anisotropic 
source. This general form includes the known cases of col- 
lisionless radiation, (4D) gravitational waves, electric and 
magnetic fields, strings and walls We assume that the 
KK anisotropic stress on the brane behaves qualitatively 
in a similar way, i.e.. 



(16) 



where ^Cf^yC^^^ is 0(1), while p'^'^ is perturbatively 
small relative to the matter energy density p, as shown 
by Eq. (p^). The energy conservation equation (|^) and 
the shear evolution equation (^ imply that 

f-(7-|)er-(7^''C^,r, r^p'^Vp, (17) 



s^. = -(e + QM)s^,, -f {k'p/ZH) rC^^y . 



(18) 



It is important to notice the special situation that arises 
in Eq. ( p7| ) when the perfect fluid background is radiation 
(7 = |). In this case the stability of the isotropic solution 
is determined at second order. Linearization about the 
isotropic expansion would lead to a single zero eigenvalue 
associated with the shear eigenvalue. 

Consider first the simpler case when 7 < | . The second 
term on the right-hand side of Eq. ([l7| ) can be neglected 
with respect to the first and the evolution of the KK 
energy density is 



p^^ = 7V(f)[t(t + /3)] 



-4/37 



(19) 



Although this analysis appears to hold for the models 
with 7 > |, it does not. In these cases r > 0, and the KK 
anisotropic stresses would have a gravitational effect that 
grows with time, invalidating the assumption of Eq. ( p^ ) 
that the volume expansion is well approximated by that 
of the isotropic solution. The solution for the shear is 
obtained from the solution of Eqs. (^ and (|l|): 

+ k'N{x)C,,4x) J[t{t + P)]-^/^^dtY (20) 

In the low energy regime at late times, t ^ /3, 

= S,.i-2/^ + i2:^!lArC,,i-(8-3'r)/37 . (21) 
07 — z 



This has the same form as the general relativity re- 
sult when anisotropic matter stresses are present 
When the KK anisotropic stress is absent, the solu- 
tion is determined by the rapidly falling mode that 
is familiar from studies of simple anisotropic Bianchi 
Type I universes with isotropic stresses. However, when 
KK anisotropic stresses are present, so that Cf^u ^ 0, 
these stresses slow the decay of the shear because of the 
anisotropic effect of their pressures. The shear evolu- 
tion becomes increasingly dominated by the C^v mode 
at late times. Note in particular that during a dust- 
dominated (7 = 1) era the two modes evolve as a^i, — 
S^^i"^ + Sn'^NC^^t-^/^. The KK mode dominates at 
large t whenever 7 > |. 

Next, we consider the radiation-dominated solution 
(7 = |). This is physically the most relevant but is 
mathematically distinct. The variables O and p take 
their isotropic universe values, and Eqs. (17) and ( [T^ ) 
have a solution of relaxation form, with Sni/ — > 0, 



K^pkk ^ [C„^(f)C"'3(f)t(f 



X [F{x) 4 
AK^p^^it,x) 



ln(4i2 +4/3i-/32)]-i 
t{t + f3){2t + (3) 



4^2 + 4/3t - /32 



(22) 
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At low energies, t P, 

K^p^^ = [2t^ Cc^pC^ Ht/U)] , (24) 
CT^. = C^, [tC^pC^ ln(tA,)]"' , U = U{x) , (25) 

where is some spatially-varying initial time, with 
t > t^. This has the same form as the general rela- 
tivity result for matter anisotropic stress |q|. For a di- 
agonal metric with expansion scale factors ai(t), it leads 
to ai{t) cx i^/2Qnt)"', where the constants are of or- 
der unity, such that Srii = 0, and are determined by the 
eigenvalues of the symmetric tracefree matrix Cy , which 
specifies the KK anisotropy shape. For example, with an 
axisymmetric (oi = a-i ^ 03) shear anisotropy, we have 
ai(i) (X a2(t) oc ii/2(lnt)i/4 and az{t) oc ti/2(lnt)-i/2. 

In the radiation case we note the distinctive slow evo- 
lution of the shear anisotropy and non-additive perturba- 
tion to the scale factors in the presence of KK anisotropic 
stress. By Eq. (p3|), the ratio of shear to Hubble expan- 
sion evolves as 



(26) 



\F{x) -f ln(4t2 + 4/3i - -f 4/3t - /32) 

so that at late times it falls only logarithmically in time: 



l/[81n(i/t*)] , t~>ii>t^ 



(27) 



In the absence of KK anisotropic stresses, s would fall off 
as s (X [f(t-h/3)]i/V(2i + /3) ^ t'^'"^- Even if s were 0(1) 
close to the string or Planck scales, it would have become 
observationally insignificant by the epoch of primordial 
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nucleosynthesis {t^s ~ 1 — 100 s) or last scattering of the 
CMB {tis ~ 10^3 

The inclusion of the bulk graviton anisotropic stress 
, which has been largely neglected in earlier studies of 
anisotropic braneworld evolution, completely changes the 
picture of long-term evolution. This anisotropic stress 
completely determines the evolution of small expansion 
anisotropics. Physically, the isotropic stress tends to 
isotropize the expansion, while the anisotropic stresses 
tend to resist this isotropizing effect. Since to first or- 
der they are both radiation fluids, it is the second-order 
effect of the KK anisotropy that dominates. Its loga- 
rithmic influence reflects its second-order character. In 
general relativity, this behaviour is only possible if there 
are anisotropic stresses from matter fields or (4D) grav- 
itational waves. In the braneworld, we can get this be- 
haviour even in the absence of tt^i^, because 5D graviton 
stresses imprinted on the brane can play the role of tt^^. 

The most interesting feature of the radiation-era solu- 
tion is that the slow decay of the shear anisotropy allows 
the shear distortion to have potentially observable conse- 
quences. The most sensitive effect is on the CMB rather 
than nucleosynthesis (this would not be the case if the 
shear decay was a power of time rather than logarithmic). 
If gravity in the bulk induces an initial anisotropic stress 
on the brane, s(tin), then during the radiation-dominated 
era, s oc [ln(i/t*)]^^, with t.^,(x) time-independent, by 
Eq. (|27|). During the short interval of dust-dominated 
evolution from the equal-density epoch t^q until tis, we 
have s cx i'^/^ by Eq. (|l|). Thus the CMB temperature 
anisotropy on large angular scales will be 



^-^(ils)«^(ti„) (1^^ 



^2/3 


'ln(tin/^*)" 




ln(tcqA*)_ 



(28) 



Crucially, the magnitude of AT/T is only logarithmi- 
cally dependent on ti„ and t*. Using 1 -I- zis = 1100 and 
1 -I- Zeq = 2.4 X W^UoHq, where Hq is the present to- 
tal matter density of the universe in units of the critical 
density and ho is the Hubble parameter today in units of 
100 kms~^ Mpc~^, we find 



AT 



4.6 X W'^n^^h''^ s{ti, 



ln(<in/<* 



ln(tcq/i*) 



(29) 



The time tia could reasonably be taken as the 5D Planck 
time ^5. By Eq. (||) the 5D Planck mass is subject to 
Ms > 10^ GeV, so that < lO^^U, where <4 « 10-*^ sec 
is the 4D Planck time. For ~ ^4 w lO"**^ sec and 
tin ^ lO^'^t*, the logarithmic term would be '~ i, and 
is relatively insensitive to quite large changes in these 
quantities. Thus, using Eq. (pT|), 



s(ti, 



(30) 



This is a much tighter constraint on the initial anisotropy 
than obtained from nucleosynthesis, Eq. (|l^). The ob- 
served large-angle temperature anisotropy in the CMB 



may have been contributed by bulk graviton effects in 
the very early universe if they have an initial amplitude 
of ^ lO^^rio/iQ- This anisotropy level is too low to have 
an observable effect on the output from primordial nu- 
cleosynthesis. 

We have shown that the tidal stresses induced on a 
perfect-fluid braneworld by bulk gravitons are the domi- 
nant factor in determining the evolution of its anisotropic 
distortion. During the radiation era this distortion falls 
only logarithmically in time relative to the expansion rate 
of the brane and can contribute a significant component 
to the large angular scale anisotropy of the CMB. An ini- 
tial (t/8 ratio of only s(tin) ~ 10"^flohQ would allow the 
observed CMB signal to be a relic mainly of KK graviton 
tidal effects. 
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